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AT NEARLY EQUAL STREAM VELOCITIES 
"by Robert G. Ragsdale and Oliver J. Edwards 
Lewis Research Center 
ABSTRACT 

An analytical and experimental study is made of the isothermal tur- 
bulent mixing that ensues when one gas in injected coaxially into a 
stream of another gas that is moving in the same direction at a comparable 
velocity. The problem of turbulent coaxial mixing of dissimilar gases has 
been the subject of a number of recent studies. These investigations have 
resulted in a number of proposed formulations for the dependence of tur- 
bulent viscosity on geometry, flow, and physcial-property parameters. 

The various expressions for turbulent viscosity are compared in this 
paper on a consistent basis. It is shown that the different equations 
predict essentially the same eddy viscosity within certain ranges of 
stream velocity ratios and density ratios. Considerable divergence of 
results occurs, however, when a given equation is applied beyond the 
range of experimental conditions for which it has been verified. 

Some of the proposed expressions include an axial dependence of the 
turbulent viscosity, others do not. Previously published data are com- 
pared with theory for three different assumed variations of eddy viscos- 
ity with axial position. It is shown that, for nearly equal stream veloc- 
ities, the eddy viscosity can be taken as constant in both the radial and 
axial directions. 
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The differences and similarities of the various analytical treat- 
ments to date have resulted, primarily, from attempts to include the case 
of equal stream velocities within the framework of Prandtl ' s original 
hypothesis for free turbulent flow. This phenomenological model predicts 
no turbulence if the two streams move at the same velocity. Experimental 
measurements have shown that this is not the case. This is not surprising, 
since, as has been suggested, some degree of initial turbulence will be 
present in the two streams. At nearly equal stream velocities, the effect 
of this initial turbulence should be most apparent. 

The results of an experimental study of the turbulent coaxial flow 
of a bromine jet into an air stream are presented. It is shown that the 
initial turbulence in the two streams can be of significant importance 
when compared with that induced by differences in the stream velocities. 
Photographs of the bromine stream show that, for equal stream velocities, 
the flow changes from a laminar appearance to a turbulent one as the ini- 
tial Reynolds number of the bromine stream is increased from 1840 to 3220. 
When honeycombs with l/8 inch diameter cells that are 2 inches long are 
placed in both the air stream and the bromine stream, the turbulent ap- 
pearance at the higher Reynolds number is greatly reduced. Photographs 
of the bromine stream for initial air-to-bromine velocity ratios of 0.85 
and 1. 52 indicate that the turbulence created by this velocity defect is 
much less than that initially present in the two streams. 

It is concluded that the turbulence initially present in the two 
streams plays an important part in the turbulent coaxial mixing of dis- 
similar gases at nearly equal stream velocities. It is therefore probable 


that simple, empirical modifications of Prandtl's expression for eddy- 
viscosity do not adequately describe the mixing process, and that addi- 
tional terms should he included to account for the initial turbulence 
present in the two streams. 

HTTRODUCTION 

Turbulent shear flow has remained a subject of interest in the 
field of fluid mechanics for a considerable number of years. Free tur- 
bulence most commonly occurs in jet and wake flows. Most of the ini- 
tial attention to jet flow was directed to the situation where an incom- 
pressible fluid issues into a quiescent environment of the same fluid. 
Considerable success was achieved by the application of phenomenological 
theories, notably Prandtl's mixing length hypothesis, and similarity 
solutions. J 

A more complex situation arises when the medium into which the jet 
exhausts is not at rest and is not of the same material. Such a system 
has been the subject of a number of recent studies. These investiga- 
tions have been prompted by interest in a gaseous-fueled nuclear rocket 
o 

engine, where a low-velocity fissionable gas is injected coaxially into 
a high velocity hydrogen propellant stream, and a supersonic combustor, ^ 
where high velocity hydrogen issues into a parallel stream of oxidizer that 
is flowing at a comparable velocity. Both situations involve the turbulent 
coaxial mixing of dissimilar gases. 

Up to a point, the approaches to this problem have been the same. 

The diffusion equation and the Navier-Stokes momentum equations are written 
for isothermal, axisymmetric, boundary- layer flow, along with the continuity 
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equation. This equation set is then applied to turbulent flow by assum- 
ing that the molecular transport coefficients can be replaced by or added 
to their turbulent counterparts. These equations are then solved by a 
transformation to a stream-function, axial-coordinate plane. All of the 
theoretical works have assumed that the turbulent transport coefficients 
are constant in the radial direction. An experimental study of hydrogen, 
helium, and argon jets issuing into an air stream indicates that this 
assumption is reasonably good) the eddy viscosity was found to decrease 
to 0.8 of the centerline value at the half -radius of the jet.^ 

To complete the analytical description of the flow field, it is 
necessary to make some algebraic statement as to the dependence of the 
eddy viscosity and the eddy diffusivity on pertinent geometry, flow, and 
physical -property parameters. It is in this regard that various approaches 
have been suggested. What is required is the equivalent of Prandtl's 
hypthesis, which stated that, in a region of free turbulence, the eddy 
diffusivity is proportional to the width of the mixing zone and to the 
difference between the maximum and minimum velocities across it. Two 
problems arise if this formulation is applied as stated. The first is that 
no turbulence is predicted for the case of equal stream velocities, al- 
though it has been observed in experimental studies. It has been demon- 
strated, however, that an eddy viscosity proportional to a velocity de- 
fect can be used to correlate turbulent coaxial mixing of dissimilar gases 

CT o 

if the stream velocities are not equal. 9 A proposal has "been made for 
the situation of nearly equal stream velocities to modify Prandtl's 
original hypothesis to include the difference in stream densities. One 
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suggestion is to replace the velocity difference with a mass flux dif- 

'Z 

ference;^ such a formulation has shown agreement with experimental data 
over the ranges investigated. This expression, however, is not alto- 
gether satisfactory, since it predicts no turbulence when the mass fluxes 
of the two streams are equal. An experimental study of this particular 
flow condition has shown that turbulence does exist for equal mass 
fluxes; the author of reference 7 proposes an expression to eliminate 
this anomaly in which the eddy viscosity is taken to "be proportional to 
the sum of the mass flux and the momentum flux. This expression is also 
shown to be in agreement with some experimental data. 

The second problem that arises in attempting to apply Prandtl’s 
free -turbulence expression to the coaxial mixing process results from 
the fact that it attempts to attribute all turbulence to the velocity 
difference between the two streams. The original equation, as well as 
all of the proposed modifications discussed above, requires that the 
eddy viscosity in the coaxial mixing region be proportional to some dif- 
ference between the two streams. This does not account for any uurbuience 
that is initially present in either of the two streams. It has been sug- 
gested that this "preturbulence " may be the dominant factor if the two 

O 

streams are at nearly equal velocities. The possible conoribution of 

initial stream turbulence and boundary layers has also been mentioned in 

357 

a number of the recent studies. * ’ Experimental evidence that initial 
stream turbulence can affect the coaxial mixing process is reported in 
reference 9, where it was found that honeycomb flow straighteners 
significantly reduced the preturbulence. This is in accord with studies 
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of the effect of grids on the eddy-diffusion coefficient in turbulent 
duct flow, where it has been found that grids appreciably reduce the 
scale of turbulence. ^ 

These two aspects of the turbulent coaxial mixing of dissimilar 
gases have been investigated and are discussed herein. A number of ex- 
pressions have been proposed for the eddy viscosity variation. Though 
the algebraic formulations appear to have significant differences, each 
expression has shown agreement with experimental data, at least for the 
range of data investigated in each case. These various relations for 
eddy viscosity are compared here on a consistent basis in order to dis- 
close their similarities and differences. Previously published data are 
compared with theoretical calculations in order to determine the axial 
dependence of the eddy viscosity. The theoretical calculations are made 
with a computer program^ > ^ that solves the axisymmetric boundary- layer 
equations with no similarity assumptions, and that incorporates arbitrary 
variations of eddy viscosity in the axial direction. The data are com- 
pared with the analysis for an eddy viscosity that increases, decreases, 
and is constant with axial position. 

Results are also presented of an experimental study of the effect 
of preturbulence on the coaxial mixing process at nearly equal stream 
velocities. Photographs of a bromine stream exhausting into a surround- 
ing air stream for various flow conditions are shown. The initial veloc- 
ity ratios, air to bromine, were maintained between 0.988 and 1.009 to 
minimize the contribution of velocity defect to the free turbulence. 
Bromine Reynolds numbers were varied from 1840 to 3230. These flow con- 
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ditions were repeated with l/8-inch passage diameter honeycomb sections 
2 inches thick in both the air and bromine streams at the injection point. 
In order to determine the relative contribution to turbulence of a veloc- 
ity defect, initial velocity ratios were varied from 0.85 to 1.5 at a 
constant bromine Reynolds number of 2300. 

EXPERIMENTAL APPARATUS AND PROCEDURE 

Figure 1 shows a schematic layout of the experimental apparatus, 
which consists of a vertical test section, metered air and bromine sup- 
plies, throttled vacuum exhaust, and electro-optical data instrumenta- 
tion. The central feature of the apparatus is the rectangular vacuum- 
tight test section, which is 8 inches square in cross section by approxi- 
mately 9 feet in height. Figure 2 shows the test section and some of the 
associated instrumentation. Two opposing faces of the test section con- 
sist of optical windows; through these is projected a collimated beam of 
white light that is intercepted by a photodetector on the opposite side. 

A narrow band-pass optical filter on the detector face passes monochromatic 
light at the peak of the bromine absorption band at 4150 angstroms. Th6 
light source and detector are rigidly mounted on a U-shaped frame to 
maintain alignment. The output of the detector provides a measure of the 
average bromine density along a given chord through the bromine stream. 

The bromine boiler is seen at the right of figure 3; bromine flow is 
controlled by the power input to an internal quartz-jacketed pancake heater 
immersed in the bromine tank. The test section is operated at the vapor 
pressure of bromine at room temperature, about 4.5 psia. Thus, the heater 
supplies the heat of vaporization required to attain a desired flow. The 
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monel toiler is coated inside with Teflon. Because of the extreme cor- 
rosiveness of Bromine, only glass or Teflon is in contact with the vapor 
until it reaches the top of the test section. There it enteres a 1-inch- 
diameter monel tube from which it is injected into the air stream. Both 
gases flow from top to Bottom through the test section. 

Lucite tube bundles at the top and bottom of the test section 
eliminate any large-scale flow oscillations in the air stream. Air and 
bromine flow rates are measured with rotameters. The sizes of the test 
section and the bromine tube were chosen so as to provide operation at 
Reynolds numbers below and above 2000, as well as velocity ratios, air 
to bromine from 1.0 up to about 50. 

The first step in running procedure is to set the desired air flow 
rate. This is done by use of an upstream-flow-control valve and a down- 
stream valve that throttles to a vacuum exhaust system. After the 
desired air flow is established at the vapor pressure of bromine, the 
bromine flow is initiated by supplying power to the boiler. The bromine 
response is virtually instantaneous, since little heat loss is incurred 
by operation at room temperature. Typically, steady flow conditions can 
be achieved in 5 or 10 minutes. 

THEORETICAL CONSIDERATIONS 

The basic analytical procedure used here to compute turbulent co- 
axial velocity and concentration fields is described in references 2, 5, 
and 11, and only the pertinent features will be reiterated. The equation 
set is composed of the continuity, diffusion, and momentum equations 
written for isothermal, axisymmetric, boundary- layer flow. The bouyancy 
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term is included, and no linearizing or similarity assumptions are made 
so that the results apply equally well near the jet origin. A 
von Mises transformation to a stream- function axial-position coordinate 
set is employed in the numerical solution. The ratio of eddy viscosity 
to molecular viscosity, (pc/p), is assumed constant in the radial 
direction, and is varied in the axial direction according to the arbi- 

_ p 

trary function, A + Bz , where z is the axial distance from the jet 
origin. The turbulent transport coefficients, pe and e, are added to 
their molecular counterparts, p and D^>, respectively. The eddy dif- 
fusivities for momentum and mass transport are assumed equal. The 
model of the coaxial flow field and the pertinent varables are shown in 
figure 3. 

In order to compare the various expressions for eddy viscosity, 
they must be rewritten in the same form. From a study of air-bromine 
coaxial mixing, reference 5 obtains the following equation: 

U € /U e W 2 

— = 0.0172 ( — - l] Re. - 250 (1) 

H VUj J J 


For Reynolds numbers that are large with respect to the constant 250, 
equation (1) can be written in the form 


(pe)<£ p^. 

p jVj p* 



( 2 ) 


In equation (2), (pe)^ is the centerline, or jet, value of the eddy 
viscosity. 


Reference 3 suggests the following expression for eddy viscosity 

(P E >« * 0 ' 025 5/2 (p e U e ' P*V (3) 
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This can be written at the jet origin as 


(p £ )£ P 1 


— ± = 0.025 

p .0 .r . p 
0 J 


V u j P J 


This expression has shown agreement with air-hydrogen data. For equal 
temperatures, the density ratio on the right side of equation (4) becomes 
the molecular weight ratio of hydrogen to air, 0.069. It should be noted 
here that the preceding expression produces an eddy viscosity ratio that 
varies in the axial direction, since both the centerline density and 
velocity are axial functions. For the purposes of this comparison, equa- 
tion (4) is used to evaluate the eddy viscosity at the jet origin, where 
the jet density and velocity are at their initial values. 

The expression proposed in reference 7 for the eddy viscosity is as 
follows : 


(p€) £ JL^e 

p.U.r. " °” 025 r. 1 p.U. + 2 

J J J J \ J J P.U. 

\ J J 


in a form similar to equations (2) and (4): This equation can be written 


at the jet origin 


(c>€) t P 1 / U e 

-VT--A. 0.025 1 + hr 

J J J e \ J 


where (pe)^ is the value of the eddy viscosity in the jet stream at the 
injection point. 

Equations (l), (3), and (5) have all been proposed to express the 


functional dependence of the turbulent viscosity pe. Obviously, they are 
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not of the same form; yet each has been shown to agree with experimental 
data. By rewriting the equations in the forms given by equations (2), 

( 4)^ and (6), it is possible to compare the various expressions on a con- 
sistent basis to see how similar or different they are. 

In reference 5 it is suggested that a viscosity ratio of the two 
streams should be included in the eddy viscosity. In order to evaluate 
this idea, equations (2), (4), and (6) can be rewritten by adding the 
viscosity ratio to the left side of the equations, and multiplying the 
numerical coefficients on the right side by the actual values of the 
ratios of the gases used in the experiments related to each expression. 
By using the viscosity ratios of the gases studied by each of the in- 
vestigators, equations (2), (4), and (6) can be written as follows: 


(P e ),£ 

p .U.r . 
0 10 


P. H /U 

-A — = 0. 228 [ =£ 

p |i . I u - 

e 0 \ 0 



(P«>t 

p.U.r. 
0 0 0 



— = 0.051 




(P€) £ p /e 

p jVo p e ^0 



(pe). P P 

£_ _e 

p .U .r . p li . 
P 0 0 0 



( 7 ) 


( 8 ) 


(9a) 


(9b) 


Equations (9a) and (9b) are both rewritten forms of equation (6) with vis- 
cosity ratios of hydrogen-air and CO^-air, respectively, since both of 
these systems were studied in reference 7. 
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DISCUSSION OF RESULTS 

In reference 5 it was shown that good agreement between theory and 
experimental data was obtained by assuming that the ratio of the 
turbulent -to- laminar viscosity p e/p was constant over the entire flow 
field. Since for the air-bromine system studied, the viscosity ratio 
was only 1.22, this assumption also results in an eddy viscosity pe 
that is essentially constant. In references 3 and 7, the proposed ex- 
pressions (eqs. (3) and (5)) yield an eddy viscosity that varies in the 
axial direction. To check the importance of an axial dependence of eddy 
viscosity, the data reported in reference 5 have been compared with the 
analysis of reference 11. The arbitrary variations of pe/p considered 
are shown in figure 4. The constant value of 6 is the one reported in 
reference 5 as best representing the experimental data for an initial 
air-to-bromine velocity ratio of 1.25, a bromine Reynolds number of 870, 
and an air Reynolds number of 1720. The other two variations considered 
were a turbulent-to-laminar viscosity ratio that is proportional to 
z 1 / 2 and one that is proportional to z' The coefficients shown for 

these two cases are those that best represented the data shown in figure 5. 
Figure 5(a) shows the comparison of the experimental data with theory for 
the three cases. The ordinate is the average bromine concentration 
normalized to the first data point. Figures 5(b) and 5(c) show similar 
comparisons for initial velocity ratios of 0.97 and 0.83. 

These results indicate that, although an axial variation of the 
turbulent-to-laminar viscosity ratio does fit the data, it is not necessary. 
A constant value is adequate, if not better. This is in accord with the 
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the case of a circular jet issuing into a quiescent environment of the 
same fluid; for this situation, it has teen established that the 
kinematic eddy viscosity is indeed constant over the entire flow field. 

Figure 6(a) shows a comparison of the various expressions for the 
turbulent viscosity as given by equations (2), (4), and (6). The data 
points on the curves from references 3 and 7 indicate the velocity 
ratios at which the analysis has been compared with experimental data. 
With the exception of the point at a velocity ratio of 2.8 (from ref. 3), 
the various expressions are in general agreement. This is quite re- 
markable, in view of the differences in the algebraic formulations, and 
the wide variations of the experimental conditions upon which they are 
based. The expressions of references 5 and 7 predict turbulent viscos- 
ities at the jet origin that are quite close; the data of reference 5 
were obtained with an air-bromine system and jet Reynolds numbers from 
255 to 3850, while the data of reference 7 were for a hydrogen-air and 
a COg-air system at jet Reynolds numbers of the order of 1 million. The 
limits of ±25 percent shown indicate the spread of the data of refer- 
ence 5. It is of interest to note that the agreement between the expres- 
sions of references 5 and 7 would not exist at velocity ratios beyond 
about 3.5. The correlation of reference 7 predicts turbulent viscosities 
that are considerably in excess of those measured in reference 5, if the 
equation is applied much beyond the range in which it has been experi- 
mentally verified. This is due to the contribution of the momentum flux 
term in equation (5), which contains a squared velocity term. 
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Figure 6(b) shows a similar comparison except that the viscosity 
ratio of the two streams is included, as given by equations (7), (8), 
and (9). The trend is to move the expressions closer together, but the 
effect is slight, since the viscosities of the gases involved do not 
differ greatly. 

The general conclusion suggested by figure 6 is that the modifica- 
tions of Prandtl's original formulation that have been obtained by 
introducing mass and/or momentum fluxes have resulted in expressions 
that are more different in algebraic structure than in actual numerical 
fact. 

Since considerable effort has been devoted to correlating the eddy 
viscosity in coaxial mixing, it is pertinent to inquire into how much 
of this turbulence is actually induced by differences between the two 
streams relative to that which is initially present. At nearly equal 
stream velocities the contribution of the preturbulence should be more 
readily detected. A series of test runs were made on a bromine jet ex- 
hausting into an air stream to investigate this effect. Photographs 
were taken of the bromine stream for a number of flow conditions, both 
with and without honeycomb sections in the two streams. Table 1 sum- 
marizes the conditions investigated. 

Figure 7 shows the bromine flow for an initial velocity ratio of 
1.009, a bromine Reynolds number of 1840, and an air Reynolds number of 
2130. This clearly demonstrates that at nearly equal strem velocities, 
a segregated laminar-like flow pattern exists at low Reynolds numbers. 
Figure 8 shows the flow pattern for a velocity ratio of 0.988, a bromine 
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Reynolds number of 3230, and an air Reynolds number of 3660. Here the 
nature of the flow is markedly turbulent, though the velocity ratio is 
essentially unchanged. This shows that, for these flow conditions, 
turbulence can he induced in the coaxial mixing region by increasing the 
stream Reynolds numbers at constant velocity ratio. 

To study the effect of upstream turbulence further, these flow 
conditions were repeated with honeycomb flow passages in both the air 
stream and the bromine stream at the injection point. An individual 
passage in the honeycomb was l/8 inch in diameter and 2 inches in length. 
Thus, the Reynolds number of the bromine stream was reduced by a factor 
of 8 upon entering the honeycomb, while the air Reynolds number was re- 
duced by a factor of 64. 

Figure 9 illustrates the flow patterns at the same low Reynolds 
number (prior to the honeycomb sections) conditions as figure 7. This 
shows that the presence of the honeycombs did not add any significant 
degree of turbulence, since a smooth, segregated flow was again obtained. 
Figure 10 shows the nature of the flow when the Reynolds numbers of the 
flow are increased as before. Here the flow is turbulent, but the level 
of the turbulence is, qualitatively, much less. Comparison of figures 
8 and 10, which are for identical flow conditions except for the honeycombs, 
shows that the honeycombs do significantly reduce the initial turbulence, 
though they do not eliminate it. 

To assess the contribution of a stream velocity difference to turbu- 
lence relative to that initially present, the air stream velocity was 
varied while keeping the bromine stream constant. Figure 11 again shows 
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the laminar-like flow pattern for a velocity ratio of 0.987, a bromine 
Reynolds number of 2300, and an air Reynolds number of 2600. Figure 12 
shows the flow pattern when the air Reynolds number is decreased to 2240, 
producing a velocity ratio of 0.85. There is no significant change in 
the appearance of the flow. Figure 13 illustrates the flow pattern when 
the air flow is increased to a Reynolds number of 4010 and a velocity 
ratio of 1.52. Some flow disturbances are apparent, but the turbulence 
is considerably less severe than that present at a higher Reynolds num- 
ber and a velocity ratio of 0.988. 

These flow studies indicate that initial turbulence plays an im- 
portant part in the nature of coaxial mixing of dissimilar gases, and 
at nearly equal stream velocities can dominate the situation. It is 
therefore unlikely that expressions which contain only differences of 
stream parameters will meet with general success and that additional 
terms will be required to account for the additional sources of turbulence 
present in the two streams. 

CONCLUSIONS 

A comparison has been made of various suggested expressions for the 
eddy viscosity in a turbulent coaxial flow system. Some arbitrary varia- 
tions of the axial dependence of eddy viscosity have been used to compare 
theory with experiment, and an experimental study of the effect of initial 
stream turbulence on the mixing region has been conducted. For the range 
of conditions investigated, the following conclusions are indicated: 

1. An axial variation of eddy viscosity does not improve the agree- 
ment of theory with experimental data that is obtained with a constant 


value . 
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2. Modifications of Prandtl's hypothesis for turbulent shear flow 
that introduce mass and or momentum fluxes rather than velocities pro- 
duce expressions whose differences are more apparent than real. These 
various expressions predict essentially the same eddy viscosity when com- 
pared on a consistent basis, as long as they are only applied within the 
range of conditions for which they have been experimentally verified. 

3. The initial turbulence present in the two streams contributes 
significantly to the coaxial mixing process, and can dominate the situa- 
tion for nearly equal stream velocities. The presence of honeycomb 
sections Im mediately upstream of the injection point can reduce the 
turbulent mixing induced by this preturbulence. 

NOMENCLATURE 

constants 

width of mixing region 

normalized average bromine concentration 
binary diffusion coefficient 
Reynolds number, 2rUp/|j. 
radial coordinate 
half -radius 

axial velocity component 
axial coordinate 

dimensionless axial distance, z/r. 
eddy diffusivity 
eddy viscosity 


A,B,C 

b 

C* 

D 12 

Re 

r 

r i/2 

U 


pe 
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H viscosity 
p dens ity 
Subscripts : 

■£. centerline 
e external 

J Jet 
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TABLE I. - FLOW CONDITIONS FOR PRETURBULENCE STUDY 


ratio, 
U /U. 
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10 E27 
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0.988 

1.009 

0.988 
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External 

Honey- 

Reynolds 

number. 
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2130 

No 

3660 

No 

2130 
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3660 
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2600 

No 

2240 

No 

4010 

No 
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2.896 
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2.896 
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Figure 1. - Schematic drawing of air-bromine system 


Figure 2. - Experimental apparatus. 
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Figure 4. - Variations of p used in data comparison for lyUj * 1.25* Rej = 870; Re e = 1720. 










